f(R) Gravity: Prom the Pioneer Anomaly to the Cosmic Acceleration 
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We use metric formalism in f(R) modified gravity to study the dynamics of various systems from 
the solar system to the cosmological scale, we assume an ansatz for the derivative of action as 
a function of distance and describe the Pioneer anomaly and the flat rotation curve of the spiral 
galaxies. Having the asymptotic behavior of action, we propose the action of f(R) = (R + A)(l + 
ln(R/R c )/(R/Ro + 2/q)) where in galactic and solar system scales it can recover our desired form. 
The vacuum solution of this action also results in a positive late time acceleration for the universe. 
We fix the parameters of this model, comparing with the Pioneer anomaly, rotation curve of spiral 
galaxies and Supernova Type la gold sample data. 
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I. INTRODUCTION 

Recent observations of the Supernova Type la (SNIa) 
and Cosmic Microwave Background (CMB) radiation in- 
dicate that universe is under positive accelerating expan- 
sion [1]. This accelerating expansion is one of the im- 
portant puzzles of the contemporary physics. A nonzero 
vacuum energy can derive universe to accelerate however 
one can ask why it is non-zero and why it is so small [2] . 
Adding a simple cosmological constant term to the Ein- 
stein equations can also accelerate universe. However, 
the problem with the cosmological constant is that why 
the energy density of matter and cosmological constant 
are in the same order at the present time? 

A varying dark energy model can partially solve this 
problem in which the density of dark energy traces the 
density of matter from the early universe to the present 
time. Modified gravity can also provide an effective time 
varying equation of state. In these models the Einstein- 
Hilbcrt action is replaced with a generic form of f(R) 
gravity [3]. In addition to the late time cosmic expansion, 
early inflationary era also can be achieved by an extra 
term to the action, as adding a cubic term to l/R gravity 
model [4] . Modifying action not only affects the dynamics 
of universe, it can also alter the dynamics at the galactic 
or solar system scales. 

There are two main approaches of metric and Pala- 
tini formalism to extract the field equations from the ac- 
tion. Considering the non Levi-Civita connection asso- 
ciating to the manifold, we can take the connection and 
the metric as the independence geometrical quantities. 
Varying the action with respect to these two parameters 
(so-called Patatini formalism) results in the field equa- 
tions [5,6]. On the other hand in the metric formalism 
the connection is the Levi-Civita connection and we do 
variation of action with respect to the metric to derive 
modified gravity field equations. The advantage of the 
Palatini formalism is that the field equations are second 



order differential equations similar to the other parts of 
the physics. 

One of the interesting issues in f(R) gravity is study- 
ing the spherically symmetry solutions. In the case of 
Palatini formalism the solution is Schwarzschild-de'sitter 
metric with an effective cosmological constant. However 
in the metric formalism the solution of non Einstein- 
Hilbert action suffers from a low-mass equivalent scalar 
field that is incompatible with Solar System tests of gen- 
eral relativity, as long as the scalar field propagates over 
Solar System scales [7,8]. One of the solutions to evade 
the solar system tests is using action in such a way that 
reduces to the Einstcin-Hilbcrt action in the low curva- 
ture regime at the solar system and at the cosmological 
regime acts as an effective cosmological constant [9,10]. 
The other problem in this issue is the consistency of the 
spherically symmetric solutions in f(R) gravity [11] that 
will be addressed in this paper. 

In this work, we try to extract an appropriate action for 
the modified gravity through the inverse solution. This 
method has been applied in the previous works both in 
the galactic [12] and cosmological scales [13]. Here we ex- 
tend the previous works to the solar system scale, study- 
ing anomalous in the Pioneer acceleration and obtain the 
appropriate action in the solar system scale. On the other 
hand following the method proposed by Capozziello et al. 
[14] we extract an appropriate action to provide a flat ro- 
tation curve in the spiral galaxies. Both the solar system 
and the galactic scale solutions are consistent with the 
modified gravity field equations in the first order of ap- 
proximation. Finally we propose a generic function for 
the action to cover all the mentioned scales and in ad- 
dition provides a late time acceleration for the universe. 
At the end we use the observational data of the Pioneer 
anomalous in solar system, rotation curve of the galaxies 
and Supernova Type la in the cosmological scales to put 
constrain the parameters of the model. 

The organization of the paper is as follows: In Sec. 
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II we introduce the modified gravity in metric formal- 
ism. Using an ansatz for the derivative of the action, 
we solve the field equation for the spherically symmet- 
ric metric and derive the dynamics in solar system and 
galactic scales. We use the observational data in the so- 
lar system as well as the flat rotation curve of the spiral 
galaxies to constrain the parameters of the model. In 
Sec. Ill we propose a generic action where in the small 
scales it reduces to the appropriate actions in the galactic 
and solar system scales. SNIa gold sample data provides 
a compatible results with the other observations. The 
conclusion is presented in Sec. IV. 



II. SPHERICALLY SYMMETRIC SPACE 

Let us take an action for the gravity depends only 
on the Ricci scalar as f(R) where in the simple case 
of f(R) = R, it is so-called the Einstein-Hilbert action. 
For a generic f(R), there are two main approaches to 
extract the field equations. The first one is so-called 
"metric formalism" in which the variation of action is 
performed with respect to the metric. In the second ap- 
proach, "Palatini formalism", the connection and metric 
are considered independent of each other and we do vari- 
ation for those two parameters independently. In this 
work we will follow the metric formalism. 

A generic form of the action depending on the Ricci 
scalar can be written as follows: 



3=1- 
2k 



^gf{R) + S r , 



(1) 



Varying the action with respect to the metric results in 
the field equations as: 

F(R)R^ - ^f(R)g^ - (V M V„ - g^a)F(R) = kT^, 

(2) 

where F = df/dR and □ = V Q V Q . From equation (2), 
we take the trace and obtain the action in terms of F 
and Ricci scalar as: 



f(R) = ±(3UF + FR-kT). 



(3) 



By taking derivative from Eq.(3) with respect to r (ra- 
dial coordinate of the metric) we rewrite this equation in 
terms of F and R as follows: 



RF' - FR' + i{UF)' = kT', 



(4) 



where ' = d/dr and f'(R) = F(R)R'. 

Replacing f(R) in favor of F(R), we obtain the field 
equation in terms of F(R): 



1 K 



+^(V M V,F-^QF). (5) 

Following the method introduced in [15], we solve the 
time-independent spherical symmetric field equation in 
the vacuum. Let us take a generic spherically symmetric 
metric as: 

ds 2 = -B(r) dt 2 + A{r) dr 2 + r 2 dd 2 + r 2 sin 2 6d(f> 2 . (6) 

Since the metric depends only on r, one can view Eq. (5) 
as a set of differential equations for F(r), B(r) and A{r). 
For the spherically symmetric space both sides of Eq. (5) 
arc diagonal and we have two independent equation. We 
rewrite Eq. (5) as: 



FR W - V M V„F - kT^ 

K M = ■ 

9nn 



(7) 



where Kua is an index independent parameter and Kua — 
Kt u ] = for all [i and v. For the vacuum space T^ u = 0, 
Kut, — K[ r i = results in: 



2F^+rF'^-2rF" = 0, 



(8) 



where X(r) = B(r)A(r). For K[ t ] - K) 
B" 



0: 



F' IX' , 2 F' IX' 2 2 
T 2Y )B - ~r ~F "2Y )B ^ B + ^ X 



(9) 



In the case of Einstein-Hilbert action (F = 1) equation 
(8) reduces to X = 1 and equation (9) reduces to the 
Schwarzschild solution. We note that for this case equa- 
tion (4) also reduces to = identity. In generic case 
having a F function of distance or as a function 

of Ricci scalar, we can obtain the metric elements from 
equations (8) and (9). 

Here we take an ansatz of F(r) = (1 + r/d)~ a for 
the derivative of action as a function of distance from 
the center, where a is a small dimensionless constant 
(q <C 1) and d is a characteristic length scale in the 
order of galactic size. Similar to the case of F = 1 we 
use equations (8) and (9) to derive X and A. We start 
with the Eq. (8), the solution results in: 

X(r)^X (l + L)- 2 ^Hl + 2 ^) 4 -^ 1 , (10) 

where X is constant of integration and for a = we 
recover Schwarzschild metric, implies Xq = 1 

In what follows we obtain metric element B(r) by solv- 
ing the differential equation of (9) for the solar system 
scales (r <C d) and galactic scales (r > d) . Once we 
derive the metric, the Ricci scalar and the corresponding 
action can be obtained. Finally we will use equation (4) 
to check the consistency of the solution. 
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A. Solar system scale (r <C d) 

In 1998 Anderson ct al. [16] reported an unmod- 
eled constant acceleration towards the Sun of about 
ap = 8.5 x 10 -10 m/s 2 for the spacecrafts Pioneer 10 
(launched 2 March 1972), Pioneer 11 (launched 4 De- 
cember 1973), Galileo (launched 18 October 1989) and 
Ulysses (launched 6 October 1990). In a subsequent re- 
port [17] they discussed in detail many suggested expla- 
nations for the effect and gave the value ap = (8.74 ± 
1.33) x 10~ 10 77i/s 2 directed towards the Sun. The data 
covered many years staring in 1980 when due to the large 
distance (20 AU) of Pioneer 10 from the Sun the solar 
radiation pressure became sufficiently small. The data 
was collected up to 1990 for Pioneer 11 (30 AU) and up 
to 1998 (70 AU) for Pioneer 10. In this section our aim is 
to explain this extra acceleration by the modified gravity 
model. 

We assume the range of r <C d in our concern and 
neglect all the higher terms of r/d and ar/d. F(r) in 
this regime reduces to: 



F(r) = 1 - -r. 



(11) 



This expression is similar to adding the first order of the 
perturbation of the action around the Einstcin-Hilbcrt 
action. From equation (10), expanding X up to the first 
order results in X = Xq = 1. 

Using (11) in the differential equation (9) we obtain 
B(r) as follows: 



B(r) = [1 



a , 3 
d r + { 2 



In 



r d 2 



r 1 a a oc 2l \Ot 1 1 
" C % + 2d + ^ r+ # r ln h-r\ 



(12) 



where c\ and C2 are the constants of the integration. For 
the case of a = (Einstcin-Hilbert action) we use the 
Schwarzschild metric as the zero order which implies c\ = 
and C2 = —6m. Using equation (12) we can obtain the 
Ricci scalar of this metric. We keep up to the first order 
of perturbation in metric as: 



B(r) = 1 



2m 

r 



a 

d r > 



(13) 



The Ricci scalar in the spherically symmetric space for a 
generic case of X is: 



R 



1 

X 



B" + -B' + A.B- —(-B' + -B) 
r r 2 X y 2 r ' 



(14) 



where substituting the metric elements, the correspond- 
ing Ricci scalar up to the first order obtain as: 



(15) 



Now to check the consistency condition, we apply the 
metric element as well as the Ricci scalar in equation 
(4). Here OF up to the first order reduce to: 



UF=—(F" 
X K 

2a 
rd 



2 , IX' B' _.. 

-F F H F ). 

r 2 X B ' 



(16) 
(17) 



Doing Simple algebra shows the consistency of the equa- 
tion for the trace equation up to the first order of per- 
turbation. 

Now we replace r in favor of R(r) in equation (11) and 
obtain F(R) in terms of Ricci scalar as: 



n^) = i+^ 



6a 2 
Rd 2 ' 



Finally integrating (18) yields action as follows: 



f(R) = R + R ln 



R_ 

Rr 



(18) 



(19) 



where R = 6a 2 jd 2 and R c is the constant of integration. 

The equation of motion for a test particle from the 
metric can be obtained. Using the week field regime, we 
define an effective potential as: 



m a 

4>n = h — r, 

r 2d 



(20) 



where the acceleration of the particles from this potential 
is 



(V 

2d' 



(21) 



The first term at the right hand side of this equation 
is the standard Newtonian gravity, however the sec- 
ond term is a constant acceleration, independent of the 
mass. We may correspond this extra term to the Pio- 
neer anomalous and constrain it with the observed value 
of a P = (8.74 ± 1.33) x 10" 10 m/s 2 which results in 
a/d ■ 



B. Galactic scale (r > d) 

Recently some of the authors have been tried to explain 
the dynamics of galaxies by the modified gravity instead 
of assuming a dark matter halo for the galaxy [14]. We 
follow the same method to extract the action with the 
ansatz of F(r) = (1 + r/d)~ a proposed in this work. We 
assume r to be larger than the characteristic length scale 
of the model d and write F(r) = (1 + r/d)~ a for a <C 1 
as follows: 



F{r) ~ (r/d)~ a ~ 1 - a\n(r/d). 



(22) 



This action can be considered as perturbation around the 
Einstein-Hilbert action. From Eq.(10), 



3 



(23) 



We follow the same procedure as we did in the case of 
solar system to extract the metric. Using equations (22) 
and (23) in (9) we obtain B(r) as: 



B(r) = ( 2 ) a 



1 



-(l-a/2) 



e 2 r 



2(1-q/2) 



(!> \l + fa-V-"/*) + e ' 2 r 2 ^- a / 2 A , (24) 



1 - a y d' 

where e^'s are the constants of integration and e[ = 
ej(l — a) for i = 1,2. For a = equations (22) 
and (23) reduce to F = 1 and X = 1 and we expect 
to recover Schwarzschild-de Sitter metric which yields 
e[ = —2m and e 2 = yjA, where A is the cosmological 
constant. For generic case when a/0, from the dimen- 
sional analysis, the constants of the integration obtain as 
e{ = -(2m) 1 -"/ 2 and e' 2 = (A/12) 1 ~ Q / 2 . We rewrite the 
metric elements after fixing e'j's: 



1 



B(r) = 



A(r) = (1 - a) 



l-(— ) 

r 



2m U-a/2 j_ f^_l(l-a/2) 



r 



j 2„ Ar 
r 12 



(l-a/2) 



(25) 



From the metric elements we get the following Ricci 
scalar 



R(r) = - 



1 



(1 — a)r' j 



■[3c 



At 2 

(12-3a)(^) 1 -/ 2 



a { 2jn /2 

2 r 



)]• (26) 



We keep Ricci scalar up to the first order term in a and 
A, then (26) reduces to 



R ( r ) = 2 



A. 



(27) 



Again to check the consistency of the solution in this 
space, we substitute (22) and (27) in equation (4). The 
solution of this space satisfies the trace equation up to the 
first order of perturbation in terms of a. We note that 
in general, for non-perturbed case the solutions might be 
inconsistent. Here the solution is valid only up to the 
first order of perturbation. 

Eliminating r in favor of R from equation (27) and 
using (22), the derivative of action obtain as: 



F(R) = ( — \R + A\r/ 2 , 
6a 



then 



f(R) 



1 



.(^L)«/2| i? + A |l+a/2 



1 + a/2 v 3a' 
For simplicity let us write action as: 
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FIG. 1. Comparing the rotation curve of galaxy in Newto- 
nian gravity, t) a 1/r (dashed- line) with the rotation curve 
from the modified gravity in the weak filed regime (solid- 
line), see equation (32). The parameters of galaxy is taken as 
M = 10 11 M Q (mass of galaxy), a ~ 1CT 6 and d ~ 10 kpc. 



f(R) = f \R + A\ 



l+ct/2 



(30) 



The dynamics of a test particle around this metric fol- 
lows the geodesic equation in weak field regime, 



0, 



(31) 



where substituting the corresponding metric elements we 
get the following velocity for a particle rotating around 
the center of galaxy: 



— (-) Q/2 



r 



1/2 



(32) 



where we ignored Ar 2 term as it is five order of magnitude 
smaller than 2m/ r. For a = we recover the standard 
Newtonian law for the rotation velocity of a test parti- 
cle, in which m = GM/c 2 and M is the mass of galaxy. 
The extra term in Eq. (32) may provide contribution to 
the flat rotation curve. Figure (1) compares the rotation 
curve of a test particle around the center of galaxy with 
an arbitrary units in the modified and standard Newto- 
nian gravity. Here we model the mass of galaxy spher- 
ically distributed up to 3.3 kpc and obtain the rotation 
curve up to 66 kpc. For a typical spiral galaxy with the 
mass of M = 10 u Mq and at the large distances (e.g. 
r > d) from the center, v ~ 200 km s _1 which roughly 
constrain a ~ 10 -6 . In the previous section we had an 
estimation for a/d ~ 10~ 26 m^ 1 which provides the char- 
acteristic length scale of the model, d ~ 10 kpc. We note 
that while equation (32) provides a flat rotation curve 
for the galaxy but it doesn't support the Tully-Fisher 
relation. 



(29) 
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III. PROPOSING A GENERIC ACTION 



In the previous section we obtained the asymptotic be- 
havior of an action in the galactic and solar system scales. 
Those two actions could describe the observations in the 
corresponding length scales without need to a dark mat- 
ter. The action for the small scale vary with a logarith- 
mic function and for the galactic scales with a power-law 
function. This asymptotic behavior of the actions guides 
us to guess a generic action which can cover also those 
two scales. Here we propose the action of 



46- 



f(R) = R + A 



R + A 



In 



R + A 



R/R + 2/a R, 



(33) 



For the range of R ^> A and R/Rq ^ 2/a, action reduces 
to 



f(R)=R + R a \nA. 

tin 



(34) 



Comparing with (19) provides Rq — 6a 2 /d 2 . Using 
|i?(r)| = Qa/rd and R/Rq 3> 1/a, satisfies the solar sys- 
tem range of r <C d. 

On the other hand for a < 1 and R ~ Rq ~ A action 
(33) can be written as: 



f(R) = (R + A) 



a , ,iZ + A, 
2 H -Rr' 



where for small a, we write action as: 



f(R) 



(R + A) 



l+a/2 



R 



a/2 



(35) 



(36) 



For a « 1, the action reduces to f(R) = R + A. We 
expect the best parameters of the model from SNIa and 
CMB experiments should be around the ACDM model. 
To see the consistency of the this action with the matter 
dominant epoch, we let R ^ A and R^> Rq. In this case 
the action reduces to f(R) — ► R (i.e. Einstein Hilbcrt 
action) and the scale factor changes as a a t 2 / 3, with 
time. 

In what follows we put constrain on the parameters of 
the model in (35). The generic FRW equation in modified 
gravity is 



3HF + 3H 2 F - -(/ - RF) = np n 



(37) 



We use SNIa gold sample with the prior of flat universe to 
constrain TO , Qa and a. Using action of (35), equation 
(37) is written as follows: 



H 2 



A + « {jff2[ _R 

6 2 X l R + 



A 



In 



R + A 

R c 1 



R 2 



R + A 



R + 2A 
(R + A) 2 



HR\ — Hq£I t 



(38) 



=L 40 




FIG. 2. Distance modulus of the SNIa new Gold sam- 
ple in terms of redshift. The solid line shows the best 
fit values with the corresponding parameters of h = 0.64, 
Sl m = 0.31, S1a = 0.69 and a <C 10 -3 with the corresponding 

X 2 min/Nd.o.f = 1.14. 



where 3-ffo^m : 
FRW equation. 

H 



For a = we recover standard 



n A ) 



1/2 



(39) 



On the other hand variation of action we respect to the 
metric guarantee the conservation of energy momentum 
and for the matter we can write p = poa~ 3 . 

From the constrain of rotation curve of the spiral galax- 
ies in the previous section, a ~ 10 -6 , we assume a <C 1 
and this term is considered as a perturbation parame- 
ter in equation (38). We solve equation (38) by per- 
turbing the Hubble parameter around ACDM solution, 
H = + aH^\ in which obtain from equation 
(39) and H^ 1 ' is calculated from 



H (i) 



1 



4 R(°) 



■{hWrW* 



■A 

+ #(0)^(0) 



fl(°> + 2A 
" R(°)+A 

R(0) 

A) 3 In 



A 



R c 



R(°r 

+ 6H(°) 
}, (40) 



where R-M and are the zero order terms obtain from 
H(°' and H(°>. The relevant parameter for comparing 
the theoretical model with SNIa data is the luminosity 
distance Dl = DL(z;Sl m ,a,R c ,A,h) and is related to 
the distance modulus of the supernovas as follows: 



H = m — M = 5 log 10 
D L = c(l + z) 



10 pc 



dz 



o H(z;n m ,a,R c ,A,h)' 



(41) 
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where the K-correction is included in the distance mod- 
ulus of the supernovas. 

We do likelihood analysis letting the Hubble parameter 
h, the cosmological parameters f2 m and Q\, and a as the 
free parameters to find the best values. The comparison 
between the observed and theoretical distance modulus 
is done by x 2 fitting as follows: 

x 2 =^2 [ /iofcs ( z ') Z Vth(zi;tt m ,a,R c ,A,h)] 2 ^ 
i ° l 

The best value for % 2 , normalized to the number of degree 
of freedom is x 2 /Nd.o.f = 1-14 . The corresponding best 
values for the parameters of the model is: TO = 0.31, 
f?A = 0.69, h = 0.64 and a <C 10~ 3 . The constrain on a 
is consistent with the results from the rotational velocity 
of spiral galaxies, a ~ 10 -6 . Finally we should point out 
that R c is not sensitive to the SNIa data. 

IV. SUMMARY AND DISCUSSION 

In this work we tried to explain the anomalous in the 
acceleration of the Pioneer spacecraft and flat rotation 
curve of spiral galaxies in the framework of the modifi- 
cation of the gravity. We started by assuming an ansatz 
for the derivative of action in terms of distance from the 
center and did the inverse procedure to derive the met- 
ric and action of the space. In the solar system scale we 
extract a logarithmic extra term to the Einstcin-Hilbcrt 
action and in the galactic scale we follow the same pro- 
cedure and found a power law action. The solution in 
both two regimes obtained as perturbation around the 
Einstein-Hilbert action and we showed that within this 
approximation the solutions are consistent with the mod- 
ified gravity equations. We note that in generic case we 
may not find a consistent solution in the spherically space 
[11]. Finally we proposed a generic action where in the 
asymptotic regimes reduce to our desired metric and ac- 
tions in the solar system and galactic scales. For the 
cosmological scales this action provides a late time ac- 
celeration for the universe. Finally we used the pioneer 
data, flat rotation curve of galaxies and the CMB and 
Supernova Type la gold sample to put constrain on the 
parameters of the model. 
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